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INTRODUCTION 


LI The concept of competing risks (CR) was first mentioned in Bernoulli's 1760 
research on predicting the death rate in a population, see Klein (2010). 


A competing risk (CR) arises when a subject is exposed to more than one cause of 
failure at the same time, where these "causes" are competing for the failure of the 
experimental unit. 


LI Consider a life time experiment with n EWN identical units, where its lifetimes are 
described by identically distributed random variables X,,X>,...,X,. Without loss of 
generality; assume that there are only two causes of failure, one of interest and the second 
of all other failure causes, (see Klein et al. 2014). 


Xj = min(X4;, X2;) vy 1, z any ne (1) 


QO The observed sample is 
(X,, 61), (Xp, 62), eons (Xn, On), 


where 


(4, if Xy< XG C 
=|) iP ke Shoe =e 


is an indicator variable denoting the cause of the failure of the i*” individual. 


wy N, 


CENSORED S#UMuPILES 


Type-I Censoring: 


Xi:m:n %2:m:n *3:mn %X4:m:n e ) e ann T 


The joint density function of the type-I censored data: 


i-Fa@)| “| [rem 


f(%41,X2, Xn) = (n— —— 


Type-ll Censoring: 


Xt-m:n %2:mn *3:mn *4mn @ e e Sonne 


The joint density function of the type-II censored data: 


(1 — Find)" a f (Ximin) 


f(%4,X2, Xn) = (n— Cee 


Withdraw | | Withdraw { 


Progressive Type-I Censoring: 


Xin = F1 X2mm X32: T2 T3 X4:mn 


The joint density function of the progressive type-I censored data: 


m k* 
f imn X2:m:n*"" > Xm:n:m) =C | | FGimn) | |i ~ F(t;)]* 
i=1 i=1 


Progressive Type-ll Censoring: 


Xy Tut a X3:mn 


The joint density function of the progressive type-II censored data: 


f 1m X2:m:n*""> Ana) =C | | f (Xim-n) [1 a F (Xin) \* 
i=1 


X, MIMI 
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Hybrid censoring scheme 


T* = min(Xm-m-nTt) This T* = max(Xn-m-nTt), This 
is called type-I hybrid is called type-I hybrid 
censoring scheme censoring scheme. 


Type-l progressive hybrid Censoring: 


Casel: Xmmn<T | 


Experiment is Terminated at T* = min{X,-n,T}- 


Case 1: {Xi-m:n s X2:m:in — Xmimn } if Xmm:n <T 
The joint density function of the type-I progressive hybrid censored data: 


m 
f t2mw X2:m:n*"" oma =C I Qin |t — F (Xi-m:n) |* 
i=1 


Case II: T < Xmmn 


X d-1:min Xx dimin 


Case 2: {X1:mn me X2:mn ae Xa:min if Xm:min >T 
The joint density function of the type-I progressive hybrid censored data: 
d 


f mn X2:m:n)*"" »Xm:m:n) =C f iman)[1 _ F (Xi-m-n)|* [1 = F(t) ]Ra 
i=1 


Generalized Progressive Hybrid censoring scheme(GPHCS) 


This generalized PHCS modifies PHCS by allowing the experiment to continue beyond time T if very few failures had been 
observed up to time t . Under this scheme, the experimenter would ideally like to observe m failures but is willing to accept a 
bare minimum of k failures; (k < m) The terminated time T* = max{ Xx.m-n,Min{Xm-m-n T }}, see Cho (2015). 


| Case l:t <= Xpaen — Xm:mn | 


Rz 


X4<m2n X2omnmn 


Case Ile Xp nh a Aen | 


X42m:n X2:mmn 


€ase Ill: Xynmen = Xmemin — F 


Generalized Progressive Hybrid censoring scheme(GPHCS) 


The likelihood function for the generalized progressive hybrid censoring scheme may 
be written as see (Zhu 2019) 


L(x|6) 


MWik1 tT Cinen) [1 7” F (Xi-m:n) |"; if T< Xk:m:n < Xm:m:n 


. _q_yd ; 
Ms fF Ximn) (1 — F Cimn) | [1 — Fy 218 ; 
if Xk:-m:n <Tt< Xm:mn 


where k; = R; incasel, and R,,, = R7,in case Ill 
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Bayesian estimation 
modified Weibull distribution Progressive Type-lIl censoring 
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ML & Bayesian method 
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INVERTED TOPP-LEONE 
COMPETING RISKS BASED 
ON PROGRESSIVE TYPE-II 


CENSORING 


Statistical properties Inverted Topp-Leone Distribution (ITLD) 


An inverse form of the Topp-Leone distribution (TLD) was proposed on the domain 
(0,c°), see Hassen et al. (2020). 


© The PDF of the ITLD is given as: (2) 


©) The CDF of the ITLD 1s given as: (3) 


© The survival function of the ITLD is given as: 
(4) 


\) The hazard function of the ITLD is given as: 


(5) 


ITLD 


© Figures of probability density and hazard functions of ITLD for some values of 8. 


| figure) 


The following observations are displayed based on PTIIC in the presence 
of CR data: 


© Xiamen < X2mn < °° < Xmm-n denote the m observed failure times. 


° 0; € {1,2} denote an indicator function of the cause of the failure of the 
i€” individual. 


¢ R,, Ro, ..., Ry» refer to the number of units removed from the study at 
the failure times X4..n-n) X2-m-n «) Xm-m-n-: 


jJ2 


Fi 
i: TI (Tisxnsn )S2 (Li:m:n )| Bike (Li:m:n)S1 (2i:m:n)| 


#=1 21 


7 (6) 
[Its (tena (eal 3 


z—tT 


where 

¢ The constant C = n(n— R, —1)...(n— 75" R; — m+ 1), does not depend 
on parameters. 

© fi (%i-m-n) Is the PDF of X3;, fo(%j-m.7) is the PDF of X>;. 

. Sy ea is the SF of X4i, 82 (Cae) is the SF of Xj. 

¢ jJ,,j2 are the number of failures due to the first and second cause of failures, 
respectively, such that j, + j. =m. 


Estimation Parameters for ITLD in presence of CR under PTIIC 


() The log-likelihood function, denoted by“, is “a 7 “yy 
In(C) + ji In (01) + jo In (02) — (20) 4+ 202 + 1) >» In (1+ %isn:n) | 


i=l 


- min(2)+ >. In (2izmnzn) 4- (01 + O02 — 1) + In] ++ 2a3-snen ) (7) 
i=1 i=1 


wr 


7 (1 -+ 02) S [Ri (—2 In (1 i Pi-m:n) -+- In (i + 2Fi-mn))| 3 


a=1 


© Now, by differentiating the log-likelihood function (7) with respect to 8, and 65, and equating them to zero, the ML 
estimator of parameters 8, and @,are, respectively 


j2 
(1 Si K can)" fea, (8) 
1 + 2Xizmin 


ji 
(1 + | en sale 
| ai Dk erie 


and 6oyrt = 


241 ln | 


O1ML = 


Pai n| 


Estimation Parameters for ITLD in presence « of ( CR under PTIIC 


©) The log-likelihood function’s second derivatives 


with respect to 6, and @, are given by j2 aes ad 2 


ae | = - 
62) 801802 062001 


@The Fisher information matrix (FIM) of the parameters 4 ares 
6, and 6. 7 # OOF 


(©) The variance-covariance matrix of the unknown y 
parameters 6, and 6, is given by | > ~ 
0. 0 


=e ~ (11) 
bom L(Oim7 + O2mML) 
Mm 


() The two-sided 100(1 - y)% , 0 < y <1, and the ECls for the unknown parameters 6, and 62, respectively, can be acquired as 
follows: 


(12) 


Bootstrap 


Estimation Parameters for ITLD in presence of CR under PTIIC 


Estimation Parameters using Bayesian method 


() We presumed independent gamma prior va \ 
to the shape parameters 6, and @, with b%1 B22 
hyper-parameters a, > 0 andb, > 0; 7™,2(41, O2|x) = an. Aas 9S1— 1 pg2— ; e (6101 +202) (13) 
n = 1,2. Then the joint prior distribution (a1) I(a2) 
of parameters 0, and 62. _ 


© The joint posterior distribution of 


i ¥ aytji—l O41 bit) ojnt (Rit) (2 In(1+2;:men)—In(1+22;:m:n)) ) 
parameters 6, and 6, is given by (41, O9|2) Oo 0; C | 


(4) 
x. ggtiely Abt PR) 2c 


© The marginal posterior distributions | - 
of 0, and 6, are, respectively, r*(0,|z) = a 1 g Oily zag n*(0o|0) = 962 1 0202 
given by PC) ie: ,* 


where ¢y = jy + Ap, and un(x) = by + (x);(2) = 2, (Rit 1) [2n(1 + Pi-men) — In(1 + 22:-m:n)]. 


Estimation Parameters for ITLD in presence of CR under PTIIC 


eo) Bayesian estimators under SE Loss function 


(~) Bayesian estimators under LINEX Loss function | -_ h 
: v)) =—Injl+—]:h#0, and n =1,2.07) 
h Up 


(Bayesian estimators under GE Loss function 


where ¢y = jy, + a,,and un(x) = b, + (x); p(x) = we (Rj +1) [2In(1 + Fimn) — (1 + 22;-mn)]. 


Estimation Parameters for ITLD in presence of CR under PTIIC 


() The marginal posterior distributions of 6, and 


(2 
0, are, respectively, given by U2 


D'(¢2) 


q1 
Uj 


(G1) 


gg eo O22 


gst h Aim and 15(0@9|2) = 


O Let Wi = 201u1 — x3, and We = 242v2 ~ x3,,, 
then 100(1 - y)% ,0<y <1, the BCls for 6, 
and @,are obtained as follows 


where X(a-y/2,2¢,)and X(y/2,26,) respectively, represent the upper and lower y/2°" percentile points of a v2 distribution 
with degree of freedom 2¢, ,¢n = jy + Gy, and Un(x) = by + W(x); w(x) — bar (R; il. 1) [2 In(1 +. Resren) = In( A sneer |e 


Numerically Study 


A simulation study is used for estimating the unknown two parameters of ITLD in presence of CR under PTIIC. 


(_) The simulation is conducted by using Mathmatica software, N = 10,000 iteration with different sample sizes for each cause 
of failure and number of failed units are 


(n,m) = (30,20), (50, 20), (60,40), (100, 70), (130, 100). 
Assume the chosen values of the parameters of ITLD are (8; ,02) = (0.5,1). 


Assume the removal of the surviving units will be done using the subsequent censoring schemes: 
Scheme 1: R; = a if i = 1,3,..,2(n — m) — 1; where the function [x] is known as the ceiling function, 
which gives the smallest integer > x; otherwise, R;= 0. 
Scheme 2: R;= a ml if i = 2,4,...,2(n — m); otherwise, R; = 0. 


Scheme 3: R;= Oifit = 1,2,3,..,m — 1,andR, = n — m. 
Scheme 4: R,= n — mandR;= Oifi = 2,3,...,m 


From Bayesian estimates (BEs), the three loss functions, SE, LINEX(h = 0.7 and — 0.7), and GE(c = 0.5 and — 0.5), are 
used. The BEs for the gamma prior are computed according to hyper-parameters (a1, a, b;, bz), which are selected by the 
same procedures that were used by Dey et al. (2016). 


Applying Numerical Study on ITLD in presence of CR under PTIICS 


simulation result for (6, = 0.5,62 =1) ) 


The performance of the different estimators 1s evaluated in terms of MSE, Absolute biases(AB). 


gamma prior uniform prior 

fips ips fig fige P) fig fips fipe 
(h=0.7) (h=-0.7) (c=0.5) (c=-0.5) BS (h=0.7) (h=-0.7) (c=0.5) 
0.0274 0.0199 0.0125 0.0276 DLL 0.0097 0.0274 0.0435 d 
0.0467 0.0121 0.0114 0.0129 0.0116 0.0467 0.0520 
0.029 D020 0.0135 0.0285 0.0094 0109 0.0290 0.0451 7025 
0.0454 0.0109 0.0123 0.0107 0.0111 0.0454 0.0506 0.0420 
0.0235 : mO108 0.0254 0.0116 0.0082 0.0235 0.0394 0.0352 
0.0449 0.0104 0.0118 0.0104 0.0106 0.0449 0.0499 0.0423 
0.0282 0.0149 0.0302 0.0083 0.0123 0.0282 0.0443 0).0306 
0.0461 0.0114 0.0130 0.0112 0.0116 0.0461 0.0514 0.0428 
0.031 0.0144 0.0293 Jee 0.0117 0.031 0.0471 0.0277 
0).0455 0.0107 0.0122 0.0109 0).0455 0.0507 0.0420 
0.028 ' 0.0123 0.0272 0.0096 0.028 Of 0.0441 0.0309 
0.0458 0.0108 0.0122 0.0107 0.0110 0.0458 0.0510 0.0426 
0.0272 3.0135 0.0282 0.009 0.0109 0.0272 UTZ 9.0433 0.0316 
0.0460 0.0106 0.0121 0.0105 0.0108 0.0460 0.0512 0.0429 
0.0266 0.0137 0.0285 0.0089 0.0111 0.0266 0.0426 0.0319 
0.0445 0.0104 0.0119 0.0103 0.0106 0.0445 0.0496 0.0415 


(n,m) Scheme Orage fips 


v The MSE values in BEs under the gamma priors were less than the MSE values of BEs under the uniform priors 


Applying Numerical Study on ITLD in presence of CR under PTIICS 


simulation result for (6, = 0.5,62 =1) ) 


The performance of the different estimators 1s evaluated in terms of MSE, Absolute biases(AB). 


gamma prior uniform prior 
~ By 05 Oy By ~ 0» Dy Oy 
im) Siheme "ML (h = 0.7) (h = 0.7 (c = 05) (c = 0.5) MBS (h = 07 (h = -6.7 (c = 05) 
0.0537 0.0216 0.0513 0.0059 0.0261 0.0241 0.0859 0.0054 
0.0955 0.0221 0.0265 0.0219 0.0231 0.0843 0.1105 0.0857 
0.0501 ys 8.0225 0.0526 0.0067 0.0271 0207 0.0821 0.0089 
0.0925 0.0218 0.0263 0.0215 0.0228 0.0818 0.1071 0.0833 
0.0549 0.0245 0.0544 0.0088 0.0291 OE 0.0254 008 0.0041 
0.0927 0.0216 0.0261 0.0213 0.0226 0.0818 0.1074 0.0830 
0.0532 0.0251 0.0555 NON 0.0298 0.0237 0.0853 0.0059 
0.0925 0.0223 0.0270 ~ 0.0219 ) 0.0233 0.0817 0.1072 0.0830 
0.0466 36 0.0217 0.0518 0.0263 DAG 0.0173 0.0784 0.0123 
0.0927 0.0220 0.0264 0.0230 0.0822 0.1070 0.0840 
0.0555 p.0251 0.0551 0.0297 t 0259 0.0877 0.0037 
0.0939 ' 0.0219 0.0265 0.0229 0.0828 0.1089 0.0841 
0.0532 0.023 0.0524 0.0275 OE 0.0237 0.0853 0.0059 
0.0923 0.0209 0.0251 0.0218 0.0815 0.1071 0.0829 
0.0479 0.0219 0.0518 a-606 0.0265 0.0187 0.0796 0.0109 
0.0875 0.0209 0.0252 0.0218 ; 0.0775 0.1011 0.0790 


1 


2 
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v The MSE values in BEs under the gamma priors were less than the MSE values of BEs under the uniform priors 


MSE for BL(h=—0.7) of 6, estimate MSE for BL(h=—0.7) of @s estimate 


0.025 
0.020 
0.015 
0.010 
0.005 


: | : 
Schemet Scheme2 Scheme3 Scheme4 esesinias Scheme1 Scheme2 Scheme3 Scheme4 


n=30 w n=50 w n=60 0D n=100 BW n=130 BW n=150 ms n=30 wp n=50 B nm=60 0D n=100 B n=130 BD n=150 
(a) (b) 


0.012 
0.010 
0.008 
0.006 
0.004 
0.002 
0.000 


MSE for BG(c=0.5) of 01 estimate MSE for BG(c=0.5) of @2 estimate 


Schemei1 Scheme2 Scheme23 Scheme4 Schemet1 Scheme2 Schemes3 Scheme4 


mw n=30 pw n=50 Bw n=60 © n=100 B n=130 BD n=150 n=30 m8 n=50 8 n=60 0D n=100 B@ n=130 BW n=150 
(c) 


3. The BE under the GE loss function (with c = 0.5) has the minimum MSE with respect to 615¢. 
4. \In all cases, the MSE values of MLE was the worst, see Figure(4). 


MSE for different estimates of @, at n=30 MSE for different estimates of 6, at n=60 


Scheme1 Scheme2 Scheme3 Scheme4 


Schemet1 Scheme2 Scheme3 Scheme4 
= ML © BS wm BL(h=0.7) @ BL(h=-0.7) e ML © BS @ BL(h=0.7) @ BL(h=-0.7) 
@ BG(c=0.5) @ BG(c=-0.5) w@ BG(c=0.5) @ BG(c=-0.5) 


Figure (4) 


4. The MSE values for BE of 05 under the BS is the minimum, while the MSE values for BL (with h = -0.7) 
of 0,were the largest, see Figure (5). 
5. In all cases, the MSE values of MLE was the worst. 


MSE for different estimates of @2 at n=50 MSE for different estimates of 62 at n=130 


Scheme1 Scheme2 Scheme3 Scheme4 Scheme1 Scheme2 Scheme3 Scheme4 
= ML @ BS © BL(h=0.7) @ BL(h=-0.7) = ML w@ BS © BL(h=0.7) mw BL(h=-0.7) 
@ BG(c=0.5) @ BG(c=-0.5) @ BG(c=0.5) @ BG(c=-0.5) 


| Figure (5) | - 


Applying Numerical Study on ITLD in presence of CR under PTIICS 


L) The AW and CP of 6, for the ECI, Bootstrap, BCI, and HPD for gamma and uniform priors. 


(n,m) 


(30,20) 


(50,20) 


l 
2 
3 
4 
l 
2 
3 
4 
l 
2 
3 


(60,40) 


— 


¥ The AWs for the HPDs are smaller than the AWs for BCI. 


Scheme 


bootstrap 


ECI Boot - p 


AW(CP 


0.7802(0.927) 0.846(0.944 
0.790€(0.928) ) 0.842 
0.786(0.926) 8861 


SSS SS eS eS oe 


0.945) 

0.938) 
0.7905(0.932) 0.848(0.941) 
0.7875(0.920) 0.837(0.945) 
0.790.931), ) 0.835(0.945) 
0.7877(0.929) 0.844(0.948) 
0.5469(0.937) _0.569(0.946) 
0.5466(0.941) ) 0.56(0.947) 
pre oa Dhara 
0.5463(0.94) 0.573(0.948) 


941 
(0.929) 
0.7878(0.935)  0.844(0.947 
(0.937) 


Boot - t 
AW (CP 
0.8722(0.884) 
0.8757(0.886) 
0.8829(0.883) 
0.8725(0.993) 
0.8611(0.886) 
().8695(0.885) 
0.8707(0.885) 
0.8714(0.994) 
0.5741(0.914) 
0.5783(0.902) 
0.5788 (0.91) 
0.5817(0.995) 


emma prior 


BCI 
AW(CP 
0.5663(0.991 
0.5635(0.993 
0.5561(0.994 
0.5685 (0.994 
0.5598(0.993 
0.5618(0.993 
0.5578(0.994 
0.5595(0.994 
0.3886(0.993 
0.3907(0.995 
0.386(0.994) 
0.3864(0.994) 


HPD 
AW(CP 
0.5326(0.903) 

i 


015814(0.905) 
0.5326(0.901) 


0.3813(0. 


“ simulation result for (6, = 0.5,0 =1) ) > 


uniform prior 


BCI 
AW(CP 
0.7786(0.941) 
0.7804(0.942) 
0.7754(0.939) 
0.78(0.94) 
0.7777(0.943) 
0.7773(0.943) 
0.7734(0.94) 
0.773(0.944) 
0.5433(0.943 
0.5429(0.948 
0.5421(0.943 

0.5472(0.947 


HPC 


0.7419(0.904) 
0. _s - 
0.7487(0.923) 
0.7473(0.923) 
0.7513(0.928) 

) 


0.7469(0.927 
0.7475(0.931 
0.537(0.946) 
0.5337(0.937) 
0.5321 (0.93) 
0.5354(0.943) 


(0 

( 
0.7467(0.924 

( 


~ simulation result for (8, = 0.5,02 =1) 


Applying Numerical Study on ITLD in presence of CR under PTIICS 


LJ The AW and CP of 6, for the ECI, Bootstrap, BCI, and HPD for gamma and uniform priors. 


uniform prior 


gmma prior 


HPD 


Bootstrap 


ECI Boot — p Boot —t BCI BCI HPC 


(n,m) Scheme 


1.1266(0.944) 
1.1289(0.947) 
1.1274(0.945) 
1.1276(0.949) 
1.1259(0.947) 
1.1276(0.949) 
1.126(0.948) 
1.122(0.952) 
0.7777(0.947) 
0.7788(0.95) 
0.779(0.949) 


0.7774(0.952) 
¥ The AWs for the HPDs are smaller than the AWs for BCI. 


AW(CP) 
1 2371(0.944) 
1.2319(0.947) 
1.2572(0.94) 
1.2269(0.939) 
0.837(0.952) 
1.2309(0.95) 
1.2441 (0.943) 
1.2108(0.933) 
0.8069(0.954) 
0.796(0.952) 
0.7837(0.946) 
0.7139(0.941) 


AW(CP) 
1.1609(0.919) 
1.1561(0.935) 
1.1662(0.921) 
1.1513(0.931) 

1.147(0.935) 
1.1507 (0.937) 
1.1526(0.936) 

1.1509(0.94) 
0.7781(0.94) 
0.7805(0.942) 
0.759(0.938) 
0.7856(0.942) 


AW(CP) 
0.8021 (0.993) 
0.8031 (0.93) 
0.8006(0.992) 
0.807(0.993) 
0.7934(0.994 
0.8022(0.994 
0.7979(0.994 
0,8028(0.994 
0.5477(0.994 
0.5501(0.994 
0.5504(0.995 
0.5498(0.995 


AW(CP) 
0.7914(0.991) 
0.8043(0.992) 
0.7911(0.991) 
0.7972(0.992) 
0.7876(0.992) 
0.7896(0.994) 
0.7908 (0.992) 
0.7871(0.993) 
0.542(0.994) 
0.5536(0.994) 
0.5464(0.993) 
0.5425(0.994) 


AW(CP) 
1.1203(0.947) 
1.1189(0.946) 
1.121(0.947) 
1.1202(0.951) 
1.1185(0.947) 
1.1221 (0.948) 
1.1186(0.95) 
1.115(0.952) 
0.7738(0.948) 
0.7722(0.948) 
0.7735(0.95) 
0.7749(0.952) 


AW(CP) 
1.094(0.943) 
1.1028(0.945) 
1.0986(0.937) 
1.0957(0.942) 
1,0927(0.943) 
1.1006(0.946) 
1.0983(0.939) 

) 
) 
) 


0.7675(0.946 
0.7654(0.948 
0.7694(0.945 


( 
( 
| 
( 

( 
1.0946(0.942 
| 
( 

( 

0.702(0.947 


@ The actual data set was initially examined by Meintanis (2007). 


Time of failure due to Cause 1: ().0333, 0.1333, 0.2667, 0.3000, 0.3167. 
().4000, 0.4333, 0.4500, 0.4667, 0.5667, 0.6000, 0.6500, 0.6667, 0.7000, 0.8167. 
().8167, 1.1000, 1.1500, 1.2000. 


Time of failure due to Cause 2: ().0500, 0.0500, 0.1167, 0.1500, 0.1833. 
().2167, 0.2333, 0.2333, 0.2500, 0.2667, 0.3000, 0.3333, 0.4667, 0.5000, 0.6500, 


.8000, 1.0333, 1.0667. 


~ A real data sets are provided in order to examine the flexibility of the ITLD. 


~ The corresponding p-value for Kolmogorov Simrnove (KS) _ test are 0.581135 and 
0.731462 for causes 1 and 2, respectively. 


0.4 0.6 0.8 ; ; : 0.2 0.4 0.6 0.8 
Data of Cause 1 Data of Cause 2 


© From the complete data, the PTIIC sample with an effective sample size of m = 30 of outn = 37 will 
be generated. The PTIIC sample is presented in this way. 


(0.05,2),(0.05,2),(0.1167,2), (0.1333,1),(0.15,2),(0.1833,2),(0.2167,2),(0.2333,2),(0.2333,2),(0.25,2), 
(0.2667,2),(0.3,1),(0.3,2),(0.3167,1),(0.3333,2),(0.4,1),(0.4333,1),(0.4667,2),(0.4667,1),(0.5,2), 


(0.5667,1),(0.65,1),(0.65,2),(0.6667,1),(0.8,2),(0.8167,1),(1.0333,2),(1.0667,2),(1.1,1), (1.2,1). 


© This comparison is exclusively depending on different PTIIC schemes, which are defined as follows: 
Scheme 1: R, = 4, R39 = 3, and R; = 0 otherwise. 
Scheme 2: R; = poe if i = 1,3,...,2(1 — m) — 1; and R; = 0 otherwise. 
Scheme 3: 
R= {2,0,0,0,0,0,0, 2, 0,0, 0, 0,0, 0,0, 2,0, 0,0, 0,0, 0,1, 0, 0,0, 0, 0, 0, O}. 
Scheme 4: R; = Oift = 1,2,...,29; and R39 = 7. 


Bootstrap uniform prior 
Boot — p Boot —t BCI 
(2.1837,4.8662) A. ) 
2.6825 
(2.6806 6.2374) 0: ) 
3.556 : 
(2.5607,5.8058 ) 
3.2451 
(1.6515,3.6316) 
1.9801 
(1.9089,5.5698) 
(2.1191, 7-8 
5.2891 
(2.0969,6.5512) 
4.4544 
(1.5327,4.1993) 
2.6666 


© The BEs were very similar to 
the ML estimates for all 
schemes. 


(2.0914,4.2091) _(1.3913,4.4163) | (1.2747,4.2428 
(2.5509,5.2¢ (1.5836,5.2711) | (1.5836,5.2711 
2.7104 3.6875 
(2.4280,4.8834)  (2.9577,4.5346) | (1.4931,4.9701) 
te ee 
(1.7638,3.2867) _ (1.0581,3.3586) | (0.9694,3.2267) 
ee ae 
(0.6950,5.1845)  (2.3937,6.1074) | (2.2702,5.9376) 

(2.8204,7.3766) 

(2.6594,6.955) 

(1.7265,4.5156) 

_ 


4.4896 oi lot 
PAanienches schemel scheme2 scheme3 scheme4 
Ai Oo 01 O2 01 O25 Ai O2 


© The AW of the Boot-t Cl of 
0, was smaller than the AW 
of the Boot-p, vice versa 


(1.3695, 4.9388) 
3.5693 

respect to @>. 2 | (2.6997,7.3358) 
3.6609 


Para- | Sch- 
per 
1 | (1.1691,4.2161) 
a, Le 3.0470 
2 | (1.4525,5.2378) 
Le 3.7804 
4 | (0.8891,3.2064) 
belie 
1 | (2.1730,5.9048) 
3.7318 
(0.0138,6.7130) 
6.6992 
3 | (2.5455,6.9169) (0.0846,6.1145) 
4.3714 6.0299 

(1.6526,4.490) 
2.8380 


(0.5492,3.6568) 
3.1076 


(2.8204,7.3766) 
4.5561 
(3.0891, 4.6660) 
1.5768 
(1.8204, 4.6673) 
2.8469 


© The AW of the Boot-t Cls of 
0,was the largest in all 
schemes. 


ML estimate 


2.69261 


4.03892 


3.34514 


5.01771 


3.15415 


4.73122 


2.04774 


3.07161 


prior 


SE 


LINEX(h = 0.7) 
uniform | LINEX(h = —0.7) 


GE(c = 0.5) 
GE(c = —0.5) 


2.69261 
2.50097 
2.92921 
2.52495 
2.63711 


4.03892 
3.75146 
4.39381 
3.87104 
3.98322 


3.34514 
3.05585 
3.72135 
3.13684 
3.27619 


5.01771 
4.58377 
5.58202 
4.80914 
4.94851 


3.15415 
2.89529 
3.48567 
2.95774 
3.08913 


4.73122 
4.34293 
5.22851 
4.53456 
4.66597 


2.04774 
1.93438 
2.18075 
1.92023 
2.00553 


3.07161 
2.90157 
3.27112 
2.94394 
3.02925 


GENERALIZED INVERTED 
EXPONENTIAL COMPETING 


RISKS MODEL BASED ON 
GENERALIZED PROGRESSIVE 
HYBRID CENSORING 


Statistical properties GIED 


A useful two-parameter generalization of the inverted exponential distribution (IED) is called the GIED, which was 
proposed by Abouammoh and Alshingiti (2009). 


© og a) _ Ai sel 
The PDF of the GIED is given as: *e Ma (1 _p " 2,0 andi>0} © 
L 
@ The CDF of the GIED is given as: 
(6) 
@ The survival function of the GIED is given as: 
(7) 


@ The hazard function of the GIED is given as: 


ae A/* 


——$ > { Re A’ (8) 
x? (1—eA/t) ° ~ 


h(x; a, A) = 


The generalized inverted exponential Distribution(GIED) 


@ Figures of probability density and hazard functions of GIED for some values of a, A. 


B a=0.5,A=1 gy a=1A=2 : @ a=0.5,A=1 gj a=1A=2 
@ a=0.8,A=1 @ a=1.5,A =2 @ a=0.8,A =1 @ a=1.5,A =2 
@ a=0.8A=1.5 gy a=2,A=3 : @ a=0.8,A=1.5 M a=2,A=3 
@ a=1A=15 gj a=3,A=5 M a=1As1.5 ff a=3,A=5 


Figure (7) Figure (8) 


Competing Risks under Generalized Progressive Hybrid Censoring Scheme(GPHCS) | 


Correspondingly, under GPHCS and competing risks model following three types of observations are 
obtained as 


OUD pepsi O1, R), —*, do, Ro), cee CE pepe, Oi, R,x)} 
if 7%. Fi Reesds = eee (Case I) 


{(X4-m:n; O1, R1), (Xo:m:n; 02, Ro), ideas (Xa:m:n, Od; Ra)} 


if A hernv a Xm:m:n> (CaselT) 


11D tates O14, R;), oo do, Ro), cosy  Aasespeas Bes Fiza it 
if Xk:m:n < Xm:m:n <7, (Case ITT) 


where 6; = 1 when X,; < X>;(failure from Cause 1) ord; = 2 when X,; > X>;(failure from Cause 2); 
Vi = 1,2,...,m. Also, R;is the number of components withdrawn at the time of the i*” failure; R;> 0. 


¢ The likelihood function can be 


L111, 2, ie FS lien 1, AO Lemni 2, A)| 

rewritten as follows in a unified 
expression: Th [f (im:n; a, A)S(Lizm:n3 1, r)| 

ee [S (pes 1, AST Dee 92, ry) 

B 
* w* 
(str a1, A)S(T; a, we Dein Re 
where 


W “denoted total number of failures observed 
experiment that terminated time 


T* = max{ Xkem:n ) MiIN{X m:m:n Tt }}- 


k: Casel 
; Wr = g 


d: CaseIJ and B= 


0; otherwise 
first and second cause of failures, 


° j,,j2 are the number of failures due to the 


m: Caselll 
respectively, such that j; + j, = W*. 


Estimation Parameters for GIED in presence of CR under GPHCS 


@ The log-likelihood function, denoted by /*, is nivted+- minted WO — A 


oa a oe es & 


w* Ww % 
2 > bt Csi 972-55) — > In (1 —e Timm ) 
é=1 i=1 


w~ 
(a1 + aa) SCR: + 1)In (1 — ce Fimm ) 
i—1 


(ay + a2) Bln (a = or) (» _ wr 


By taking the first differentiate with respect to (@;,@,A ), and then equating them to zero. It cannot be solved directly . 


©The ACI of the unknown parameters a1, a2, and A can be obtained as follows: 


aimL~+ZyV/Var(@imz), @m~+ZyVVar(@amz), and AmiL+Z ay Var(\mz)} (0) 


) 


y 2 


and — — 


where Z,/2 is the percentile of the standard normal distribution with right-tail probability y/2. 


Estimation Parameters for GIED in presence of CR under GPHCS 


@The joint prior distributions of the 
parameters a,,a,and / are written as 
follows: 


1 1,02 1,83 
b; by bs 


fi-1 yaa—1 \as—T,—(bratnoa-Hhs4) (15) 


where the hyper-parameters a, and b,, for | = 1,2,3, respectively. 


@ The joint posterior distribution of 


parameters a, ,azand A is given by I 2 


T'(ay,09,Alt) x a 
(16) 


y ast 1 tn Lidar = Sie [ln( [Ie ema] Li mn 


where o(2,7,) = \-,(Ri +1) In(1—e-/*em") 4 B (n ~W Ley Ri) In(1-e-/"). 


we will utilize the MCMC method to obtain the Bayesian estimators. 


Numerically Study 


© A simulation study is used for estimating the unknown three parameters of GIED in presence of CR under GPHCS. 


© The simulation is conducted by using R software, N = 5,000 iteration , MCMC= 2000 sample with n-burn=500 where 
different sample sizes for each cause of failure, the number of failed units with the various minimum number of failure units 


are 
(n,m, k) = (60,40, 25), (60, 40, 30), (60, 50, 25), (60, 50, 30), (100, 70, 45), (100, 70, 60), (100, 85, 45), (100, 85, 60) 


under the different pre-fixed time T = 1.2,andT = 2. 


© Assume the chosen values of the parameters of GIED are (@1,Q@2,4 ) = (0.8,1.2,1.5). 


© Assume the removal of the surviving units will be done using the subsequent censoring schemes: 
Scheme il: R,= n — mandR; = Oift = 2,3,...,m. 
Scheme 2: R;=1 if t = 1,2,...,2 —m; otherwise, R; = 0. 
Scheme 3: R,;= Oifi = 1,2,3,..,m — 1,andR, = n — m. 


(© From Bayesian estimates (BEs), the three loss functions, SE, LINEX(h = 1.5 and — 1.5), and GE(c = 0.5 and — 0.5), are 


used. The BEs for the gamma prior are computed according to hyper-parameters 
(a4, b,) = (12. 8, 16), (a2, bz) = (28. 8, 24. 0), (a3, b3) = (45. 0,30. 0) 
(by letting the mean of the marginal prior distribution be true parameters and its variance is 0.05 for each prior) see Mohie 


El-Din et al. (2019) 


(60,40,25) 


378 


0.0007 
0.0105 


oO. 
0. 
oO. 
8 
8 


0.0011 
0.0099 


yO CO 


(60,50,25) 


; ‘ono7 | Cox 

(60,50,30) tnenee.. 6 0.0286 
0.0102 0.0124 

0.0185 0.0251 


0.0097 0.0114 


YE — 


(Simulation result for (@1,a,a) = (0.8,1.2,1. B) ) 


MSE for different estimates of a, MSE for different estimates of az 
under Scheme? at (n,r)=(60,1.2) under Scheme’ at (n,t)=(60,2) 


. m=40,k=25 m=40,k=30 m=50,k=25 m=50,k=30 m=40,k=25 m=40,k=30 m=50,k=25 m=50,k=30 


=| ML OBS @ BL(h=1.5) @ BL(h=-1.5) = ML OBS @ BL(h=1.5) @ BL(h=-1.5) 
@ BG(c=0.5) m® BG(c=-0.5) m BG(c=0.5) m BG(c=-0.5) 


| Figure (9) | 


Applying Numerical Study on GIED in presence of CR under GPHCS - 


“simulation result for(ay,a2,4) = (0.8,1.2,1.5)>) 


As t increases for fixed n, m, and k, the MSE of the different BEs increase, while the 
MSE of ML decrease, see Figure (10). 


MSE for different estimates of a, MSE for different estimates of ap 
under Schemez2 for (n,m)=(60,40) under Scheme2 for (n,m)=(100,70) 


T=1.2, k=25 T=1.2, k=30 T=2 k=25 T=2 k=30 t=1.2,k=45r=1.2,k=60 T=2,k=45 T=2,k=60 
= ML © BS m BL(h=1.5) @ BL(h=-1.5) = ML © BS m BL(h=1.5) @ BL(h=-1.5) 
= BG(c=0.5) # BG(c=-0.5) m BG(c=0.5) = BG(c=-0.5) 


| Figure (10) | 


MSE for different estimates of A MSE for different estimates of A 
under Scheme 2 and at T=1.2 under Scheme 1 and at T=2 


n=60,m=50,k=30 n=100,m=/70,k=45 n=60,m=40,k=30 n=100,m=70,k=60 
m= ML © BS m BL(h=1.5) @ BL(h=-1.5) = ML © BS @ BL(h=1.5) @ BL(h=-1.5) 
m BG(c=0.5) m BG(c=-0.5) @ BG(c=0.5) # BG(c=-0.5) 


Figure (11) 


040.40 


a ) A (i) ) h 
(087 748 0.582.672 0.7860. 69 
0.9988 0.9908 0.9962 = 0.9976 0.9906 0.944 
090 O27. 075 OPAL O55 
0,998 1.0000 0.9964 = 0.9904 1.0000 0.9950 
(609 0,739 = 0.580 = 055 0.726 (0,838 
0.9990 0.9996 0.9950 0.9872 0.9996 0.9934 
(600 O72? 0.883 O45 0.715 0.871 
0.9962 0.9096 0.9904 0.9950 0.9996 0.9950 
(604 = ,734 0.885.080 0,722.4 
09978 1.0000 0.9044 = 0.972 0.9906 0.9924 
O61 7340.49 OT 72283 
0,9986 1.0000 0.9966 = 0.9978 0.9998 0.9956 


@ The original data set was originally monitored by Xia et al.(2009), in a jute fiber-breaking strength 
experiment. In this experiment, two different gauge lengths, 10 mm and 20 mm, were used to 
influence the jute fiber’s breaking strength failure data. 


Time of failure due to Cause 1: 0.732167, 0.836, 1.68583, 1.81567, 2.051, 2.35633, 2.52467, 2.72333 
0.95417, 3.05967, 3.5355, 4.20067, 4.38167, 4.8545, 6.065, §.30717, $.88738, 6.27367, 6.3005, 7.03517 
8.14333, 8.8425, 9.84133, 10.6277, 11.1915, 11.5622, 11.679, 11.7443, 12.1205, 12.9696. 


Time of failure due to Cause 2: 0.6125, 0.759667, 0.800167, 1.191, 1.3025, 1.662, 1.8075, 1.94983, 
1.09767, 2.43967, 2.77483, 3.11883, 3.13083, 3.336, 4.0755, 4.744, 5.845, 6.2635, 6.98367, 7.61, 9.124 
9.6367, 9.69838, 9.7505, 9.90483, 11.0443, 11.4608, 11.7803, 12.6117, 12.7593 


® A real data sets are provided in order to examine the flexibility of the GIED. 


® The corresponding P-values from the Kolmogorov Simrnove (KS) test are 0.535989 
and 0.367345 for causes 1 and 2, respectively 


Fitted GIED's CDF 
Fitted GIED's CDF 


Data of Cause’ Data of Cause2 


0.732167, 0.759667, 1-191, 1.3925, 1.662, 1.68583, 1.94983, 2.051, 2.35633, 


2.02467, 2.72333, 2.95417, 3.05267, 3.11883, 3.336, 3.5355, 4.0755, 4.33167, 


4.8545, 5.065, 5.39717, 5.845, 5.88733, 6.2635, 6.27367, 6.3905, 7.03517, 
8.8425, 9.124,@-69333, 9.7595, 9.84133, 9.90483,40.6272 11.1915, 11.4693, 
11.5692. 11.679. 11.7445. 12.1905. 
Let R = (0,0,1,0,1,0,0,2,0,1,0,0,1,0,2,0,0,1,0,1,0,0,2,0,1,0,0,1,0,2,0,0,1, 0,1, 0,0, 2, 0,0), 
consider n = 60, k = 30,m = 40 and a different value of tT , then three different GPHC scheme are 
provided as below: 


LI Scheme 1: Suppose t = 6, since T < X39.49.69, then the experiment would have terminated 
at X39.49-69, With j, = 18, j2 = 12, R¢ = 0, 
R* = (0,0,1,0,1,0,0,2,0,1,0,0,1,0,2,0,0,1,0,1,0,0,2,0,0,0,0,0, 0,18). 
LI Scheme 2: Suppose t = 11, since tT > X39.49-69, then the experiment would have terminated 
at X30:40:60 < T < X40:40:60, With j,; = 20, j2 = 14, Rp = 9, 
R* = (0 ,0 ,1 ,0 ,1 ,0 0 ,2 0 ,1 ,0 ,0 1 ,O .2 0 0 1 0 1 ,0 0 2 0 1 0 0 1 0 2 0 0 1 0). 


LI Scheme 3: Suppose t = 13, since T > X49.49.69, then the experiment would have terminated 
at X 40-40-60; with It = 20; J2 = 15, R* =Rand Re = 0. 


schemes 
(| 


REAL DATA ANALYSIS scheme schemed 
| 


a fe[t [ale 


0 

Bd2 | 39492 | 0.7592 | 0.4500 | 4.6717 
208 | 0.8298 | 0.2081 | 0.1454 | 0.8803 
) 


Cl} (l} 
mi Estimate beedteadbef 
@ The BEs were very SER 01336 | 0.1010 ) 0.8016 | 0.1548 | 0.120 
ny 04246 | 0.2849 ) 3.4283 | 0.5158 | 0.3592 | 3.8597 | 0.7520 | 0.4036 ) 4.6990 
— tre 
: | 


close to the ML 
estimates in all 
schemes. 


) 

§ 

! ] 

, f 
5} | QALZL | 0.2749 | 0248 | 04995 | 0.3401 | 3.4012 | 0.7228 | 0.4304 | 4.1047 
5 18 

i 

j 


3 
| 
0108 | 0.3 
1521 | 0, | 
| 0.4999 | 0.3 
| 8802 | 0.1531 | 0.1 | 
Sirsa [INEX(h = -1.5) 13185 | 0.7853 | O60! 
standard error (SER) of SER | (130 | (0005 | 9418} 0.1833 | O10 | 0816 | 02064 | 0.1418 | 0.988 
03049 | 0.2571 | 3.2918 | 4837 | 0.3313 | 3.744 | 0.7116 | 0.4213 | 4.4748 
was the smallest. 
7} 0. Q.a00 | 0.3 
1327 0.1525 | 0.1 


i 
i 
) b 
i 
b 
BEs under uniform prior , A 
| 8 
QA147 | 0.2735 | 3.3834 | 0.5000 ) 0.3498 | 3.8221 | 0.7385 | 0.4429 | 4.5043 
0.1327 | 0.0996 | 0.7836 | 0.1525 | 0.1200 | 0.7537 | 0.2041 | 0.1413 ) 0.8362 


REAL DATA ANALYSIS | 
Para- | Sch- 
eme ACI 
1 (0.1672,0.6908) 
2 


(0.2183,0.8252) 


Bootstrap 
Boot — p Boot —t 
(0.4556, 1.8328) | (0.4541,0.7053) 
1.3772 0.2512 
(0.5634,1.6193) | (0.5598,0.8556) 
0.6069 1.0559 0.2958 
(0.7745,1.7349) | (0.7736,1.1370) 


3 | (0.3513,1.1672) 
| = | 0.8159 0.9604 0.3634 
1 | (0.0880,0.4840) (0.2878,0.5831) | (0.2877,0.4042) 


(0.3697,0.8027) | (0.3692,0.5651) 
0.4330 0.1959 
(0.4620,0.9876) | (0.4614,0.6877) 


“2 2 | (0.1284,0.6021) 
0.4737 


(0.1706,0.7404) 


(1.9344,5.0765) _ (3.5802,7.9207) 


(4.7785,7.2975) 
2.5190 


L} Cls and AW for the estimates of the real data 


Uniform prior 


BCI 
(0.2126,0.7226) 
0.5100 
(0.3451, 1.0465) 
0.7014 
(0..4215, 1.2445) 
0.8235 
(0.1286,0.5083) 
0.3797 
(0.2031,0.7165) 
0.5134 
(0.2351,0.8117) 


(3.0107,6.3617) 
3.3511 


HPD 
(0.1860,4.9962) 
0.4796 
(0.3069,0.9902) 
0.6833 


(0.1086,0.4698) 
0.3612 
(0.1865,0.6758) 
0.4893 
(0.2099,0.7378) 
0.5279 


(2.8593,6.1689) 
3.3096 


oo 
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